Unified theory of spin-dynamics in a two dimensional electron gase with arbitrary 
spin-orbit coping strength at finite temperature 
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We study the spin dynamics in the presence of impurity and electron-electron (e-e) scattering in a 
III-V semiconductor quantum well with arbitrary spin-orbit coupling (SOC) strength and symmetry 
at finite temperature. We derive the coupled spin-charge dynamic equations in the presence of 
inelastic scattering and provide a new formalism that describes the spin relaxation and dynamics in 
both the weak and the strong SOC regime in a unified way. In the weak SOC regime, as expected, 
our theory reproduces all previous zero-temperature results, most of which have focused on impurity- 
scattering induced spin-charge dynamics. In the regime where the strength of the Rashba and linear 
Dresselhaus SOC match, known as the SU(2) symmetry point, experiments have observed the spin- 
helix mode with a large spin-lifetime whose unexplained non-monotonic temperature dependence 
peaks at around 75 K. As a key test of our theory, we are able to naturally explain quantitatively 
this non-monotonic dependence and show that it arises as a competition between the Dyakonov- 
Perel mechanism, suppressed at the SU(2) point, and the Elliott- Yafet mechanism. In the strong 
SOC regime, we show that our theory directly reproduces the only previous known analytical result 
at the SU(2) symmetry point in the ballistic regime. It also explains, as we have shown previously, 
the rise of damped oscillating dynamics when the electron scattering time is larger than half of the 
spin precession time due to the SOC. Hence, we provide a unified theory of the spin-dynamics in 
two dimensional electron gases in the full phase diagram experimentally accessible. 

PACS numbers: 73.21. Fg, 72.25.Dc, 72.25.Rb, 72.10.-d 
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I. INTRODUCTION 

The spin degree freedom of electrons and the control of 
its dynamics by electric means has played an increasingly 
important role in condensed matter research and spin- 
tronics device applications over the past decades. This 
electric manipulation is achieved via the SOC, which has 
been a key factor in new spintronic devices as well as the 
emergence of new fields such as spin Hall effect^— spin 
Coulomb drag,— ~— and topological insulators. 7,8 

These fast developing SOC-based fields highlight the 
need to fully understand the basic spin-charge transport 
dynamics in the full phase diagram experimentally ac- 
cessible, both at zero and finite temperature. The spin 
relaxation and dynamics in the presence of weak SOC is 
a mature topic of research, and a theoretical perturba- 
tive treatment on the SOC gives an accurate description 
of the spin dynamics at zero temperature^— In the op- 
posite limit of strong SOC, defined as fl so T > 1 where 
fl so is the spin precession frequency due to the SOCs 
and t is the momentum scattering time, only a few the- 
oretical work o 11 ' 13 were able to derive analytically the 
spin dynamics at one specific point in the phase dia- 
gram; specifically when the Rashba SOC strength is equal 
to the linear Dresselhaus SOC strength, where a novel 
spin dynamic mode, the spin- helix mode (SHM), with 
large spin-lifetimes has been observed experimentally. e 
One theory work was able to analytically derive the spin- 
dynamics modes at zero temperature and successfully ex- 
plain the different damped oscillatory modes observed in 
this strong SOC regime^ However, until now, analytical 
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FIG. 1. Spin life time of the dominant spin-helix mode vs. 
temperature, extracted from the Ref. [TH (Green stars). The 
blue dash line is our theoretical result (see Sec. IIII C| . 

theories of the spin dynamics in 2DEGs have primarily 
focused on zero temperature . 10 ' 11 ' 13 ' 15 The rich physics 
observed in the experiments as a function of tempera- 
ture, which depict the transition from the strong to weak 
SOC and shows a peak of the spin lifetime of the SHE 
mode at ~ 75 K, has remained unexplained up to now. 15 
Here we develop a microscopic approach valid at fi- 
nite temperature and for any strength of SOC within 
the same derived expression. To show the validity of 
our method, we focus on the spin dynamics in the 
n-type III-V quantum well which has been recently 
well studied experimentally.—'^— Besides reproducing 
all known limits and previously derived results at zero 
temperature^ ' 10 ' 12 ~ — it succeeds dramatically in describ- 
ing quantitatively the non-monotonic behavior of the spin 
lifetime of the SHE mode as shown in Fig. [TJ Previ- 
ously it was speculated that this non-monotonic behavior 
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could be attributed to the increase in importance of the 
cubic Dresselhaus SOC, which decreases the SHM life- 
time. However, a simple analysis (see Sec. III. A) shows 
that in the range from zero to 100 K this term does not 
change significantly and, therefore, it fails to account for 
this behavior. Our theory shows directly that this non- 
monotonic behavior arises from the competition between 
the Dyakonov-Perel (DP) mechanism, suppressed at the 
SU(2) point, and the Elliott- Yafet (EY) mechanism. 

Our microscopic approach is based on the non- 
equilibrium Green's function formalis m 10 ' 19 which allows 
us to fully understand the temperature dependence of the 
spin relaxation in the presence of both disorder and e-e 
inelastic scattering. The e-e interaction is shown to dom- 
inate the momentum scattering time above a certain fi- 
nite temperature (35 K in the relevant experiments) and 
controls the transition between strong and weak SOC 
regime*^ The e-e scattering time in the 2DEG system 
has been theoretically discussed previously^ and, even 
though it does not affect the charge transport due to net 
momentum conservation of the system (see appendix |B|) . 
the spin-current and spin decay can be affected strongly 
by it, e.g. spin Coulomb dragj^~— By incorporating both 
spin independent scattering and SOC scattering poten- 
tial, both DP and EY mechanisms are considered in the 
frame of the quantum kinetic equation. 

Our paper is organized as follows. In Sec. |H]we intro- 
duce the quantum kinetic equation (QKE) based on the 
non-equilibrium Green's function formalism in the pres- 
ence of the three types of SOC present in III-V based 
2DEG. Within Sec. Ill Al because we want to consider 
finite temperature, we evaluate the contribution of the 
e-e interaction in the collision integral of the QKE and 
show that it dominates the momentum scattering time 
for temperatures above 35 K. In Sec. Ill B I we generalize 
the spin-charge kinetic equation from the zero temper- 
ature limit^ to the finite temperature in the presence 
of inelastic scattering. Within Sec. Ill CI we present our 
procedure which allows us to have a theory valid in both 
the weak and the strong SOC regime, with Eq. (|35|) and 
(|36f being the pivotal results that allows us to evaluate 
the full phase diagram of the spin-charge dynamics. 

In Sec. IIIII we focus on the temperature dependence 
of lifetime of the dominant SHM near the SU(2) sym- 
metry point, were we show that the thermal average 
of the SOCs strength is almost unchanged from OK to 
100K. This discards the increase of the cubic Dressel- 
haus strength as source of the non-monotonic behavior 
of the spin-lifetime as mentioned above^ In this sec- 
tion we evaluate the EY mechanism and show quantita- 
tively that the non-monotonic temperature dependence 
of the enhanced-lifetime of the SHM is the result of the 
competition of the EY and DP spin relaxation mecha- 
nism. Finally in Sec. IIVI we reproduce the spin relax- 
ation eigen-modes in the presence of the different SOCs 
which confirm that our method is also valid in the strong 
SOC regime. 



II. THE QUANTUM KINETIC EQUATION FOR 
THE 2DEG WITH THE GENERAL SOC 



The non-equilibrium spin polarization can be described 
by the non-equilibrium Green's function formalism, also 
known as Keldysh formalism. 19 This formalism was used 
in the noninteracting system with short range disorder 
and in the presence of weak Rashba SOGiS or the equal 
magnitude of Rashba and linear Dresselhaus SOCs with 
zero cubic Dresselhaus SOC.™ However, when consider- 
ing the temperature dependence of the spin life time, we 
have to generalize this method to the interacting system 
with a general SOCs. Here, we focus on the 2DEG in a 
quantum well such as GaAs/AlGaAs£±£i±L2i 

In the 2-D semiconductor heterostructures, the Hamil- 
tonian takes the form 



H=— +b(k)-<r, 



(1) 



where b(k) is the effective magnetic field and contains 
three types of SOCs, namely the linear Rashba 22,23 SOC 
and the linear and cubic Dresselhaus^ SOCs which take 
the form 

b fl (k) =a(-k y ,k x ), (2) 
b Dl (k) =Pi(k y ,k x ), (3) 
b° 3 (k) = -2/3 3 cos29(-k y ,k x ), (4) 

where hp is the Fermi wave vector. Here we take as 
the angle between the wave vector k and the [110] direc- 
tion which we take to be the local x-axis in our coordi- 
nates. The above SOCs split the spin-degenerate bands 
and dominate the spin dynamics in the 2DEG. The cor- 
responding SOC Hamiltonian takes the form: 



H so = ^\k x a v + \ik v a 



y u x 



byCJy 



b x &x, 



(5) 



where Ai = a + /3 t - 2/3 3 cos 29, A 2 = j3\ - a + 2f3 3 cos 29. 
The retarded (advanced) Green's function of the Hamil- 
tonian in Eq. ([T]) takes the form 



G R ( A \E,k) 



(E- 



(6) 



The nonequilibrium state of the system can be de- 
scribed by introducing the contour-ordering Green's 
function in the Keldysh space as 



6(1,2) 



G R (1,2) G K (1,2) 
G A (1,2) 



(7) 



where 1 and 2 stand for the condensed notation 1 = 
(xi,s z i,ti). 

In the presence of these general SOCs, the kinetic part 
of the quantum Boltzmann equation has the for m 10 ' 25 

d T G K + ~{V, -V R G K } + i[b(k) • (7, G K ] 



[(Z R G K - G K Yi A ) 



(G R Yi K - Yi K G A )\ 



(8) 
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FIG. 2. The two self energy we consider in this work. The 
double wiggle in the first diagram is the effective e-e interac- 
tion in the random phase approximation (RPA). The dashed 
line in the second diagram represents the impurity scattering. 



where V = [• • • , • • •] and {• • • , • • • } stand for com- 
mutation relation and anti-commutation relation respec- 
tively, and T, R , E A and T, K are the retarded, advanced 
and Kcldysh self energy and the corresponding Feynman 
diagrams are shown in Fig. [2] The self energy of the e-e 
interaction have the more complicated forms 26 



-A" 



(G R - G A ) o (D R — D A ) + G K o D 



K 



(9) 



where the symbol o denotes integration over all inter- 
nal energies and momenta, D R ( A ' K ^ are the full dressed 
propagator of Coulomb interaction 2 ^ and G^^ A ' K ' is the 
electron Green's function in the self energy diagram Fig. [2] 
and its energy and momentum are denoted as E' and k' 
respectively. The underline is to emphasize its difference 
to G R ( A ' K>> which is the electron Green's function out of 
the self energy; i.e. the bare self-energy. In the limit 
<i = t 2 and xi = x 2 , G K \t\ = i 2 ,xi = x 2 ) = 1 - 2n(x, t) 
where n(x,t) = ip' (x,t)i(j(x,t) is the electron density op- 
erator. 



A. Collision integral of e-e interaction 

The treatment of the collision integral contribution 
from the impurity scattering has been well established. 10 
Therefore we will focus on the e-e interaction in the colli- 
sion integral which may not be as familiar to the readers 
and which, as shown below, can dominate the electron's 
momentum relaxation time at finite temperatures, even 
if it has no direct effect on the total momentum of the 
system. The first e-e scattering term has the form 



(G K o (D 1 



D A 



{G K -G A )oD K ) G K . (10) 



-,K 



The Keldysh Green's function can be written as G K — 
G X + SG K where G${8G K ) is the Keldysh Green's func- 
tion of the equilibrium (non-equilibrium) part. The 
terms containing the first order of SG K and SG_ K in 
Eq. (fT0|) have the form 



{Go ° 



G J 



(D R - D A ) + (G R 
+6Gf o (D R — D a )Gq 

—6G K + 5G K o (D R - D a )Gq , 



D K )SG 



K 



(11) 



dE' d 2 k' 
2tt (2tt) 2 

x ( tanh( 2£r> 



{D R - D A ){G R - G A 



coth(— ) 



(12) 



with uj — E — E' . Similarly, the nonquilibrium collision 
integral in the second term on the right hand side of 
Eq. ([5]) can be written up to the first order of 6G K as 

G R 5T, K — 5T. K G A = 6G K o D K (G R — G A ) . (13) 

Substituting Eq. dTT]) , (H3) , C31) into the right hand side 
of Eq. ([5]) , the non-equilibrium e-e collision integral has 
the form 



I ee (6G K ,SG K ) = I ee (SG K )+l ee (SG K ), (14) 



where 

Iee(SG K ) 



(GfoiD 1 



D A ) + (G 1 



G J 



i — SG 



K 



D K ) 5G K 
(15) 



and 

Lee( s G K ) = SG K o [{D R - D a )Gq + D K (G R - G A )\ 



SG K o (D R - D A ){G R - G J 



E 



x r nh( 2^r j 



coth( 



2k h T 



(16) 



We note that I ee has a similar integrand as the one used 
in calculating the momentum scattering time of e-e in- 
teraction shown in Eq. [121 

Unlike the case of only considering impurity scatter- 
ing ; 10 ' 13 at finite temperature the collision integral con- 
tains inelastic e-e scatterings which scatter the electron 
to different energy states and make the collision integral 
of the right hand side of QKE Eq. ([5]) more complicated. 
Therefore we first give a brief discussion of the e-e inter- 
action in the 2DEG. In the experiments considered the 
spin splitting energy A so is much smaller than the Fermi 
energy e_F = 400 K, say A so /e_F <C lr^ and therefore the 
screening of the Coulomb interaction in this SOC system 
is treated to be the same as in the non-SOC system. The 
inverse screening length in the two dimensional system 
has the form 20 



27re 2 AT 



(17) 



where No is the density of state at the Fermi surface, e 
is the electron charge and to is the dielectric constant in 
the vacuum. For the GaAs/AlGaAs, = ^p- where 
m* = 0.065 too is the effective mass of the electron in 
the quantum well and Too is the mass of electron in the 
vacuum^ In this case, Kd ~ 3.08 x 10 s cm -1 , which is 
much larger than the Fermi wave length kp — ^Jlixno — 
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t p vs T. The blue stars are the momentum scat- 



FIG. 3. 

tcririg times extracted from of D s vs T in the supplemen- 
tary material in Ref. [l5| by using the relation r = 2D s /v% 
where Vf = — f- is the Fermi velocity and is estimated to be 
4.32 x 10 m/s. The red line is the theoretical estimate of 
T e - e based on the theory of the Ref. I2TJ . 



2.24 x 10 6 cm -1 , where n = 8 x 10 11 cm~ 2 is the den- 
sity of electrons^. Therefore the Coulomb interaction 
is strongly screened in the 2DEG considered in the ex- 
periments and we can treat the e-e interaction as angle 
independent scattering. On the other hand, the peak 
of the enhanced-life time of the SHM happens around 
75 K which is much smaller than the Fermi tempera- 
ture 400 Kf^ The non-equilibrium electrons distribute 
around the Fermi surface within the energy range of k B T 
,where kg is the Boltzmann constant and T here is the 
system temperature without confusing it with the time 
variant T in Eq. (jHJ. Therefore the non-equilibrium elec- 
trons around the turning point 75 K is in the regime 
6f ^> k B T ^> |ejt — cf\, where cf is the Fermi energy 
and £fc is the electron energy with momentum k. The e-e 
scattering time in this regime is estimated theoretically 



,.20 



~8h 



In 



eg 
k B T 



(18) 



It is noted that the e-e scattering time in Eq. (fT8[) is inde- 
pendent on the energy of the electrons and there is only 
one parameter, ep, that we need from the experimental 
data to estimate the e-e scattering time. In Fig. [3j we 
compare the momentum scattering time extracted from 
the experimental data in the supplementary material in 
Ref. [[a and the theoretically estimate based on Eq. (fT8"|) . 
They match very well when the temperature is above 
30 K. The mismatch below 30 K is because the non- 
equilibrium electrons excited by the optical field are be- 
yond the energy range of fc^T around the Fermi surface 
in the case of T < 30 K and Eq. (fT8| fails in this temper- 
ature regime. However, in a wide temperature regime, 
from 30 K to 150 K, Fig. [3] indicates that the momentum 
scattering time felt by the spin relaxation is dominated 
by the e-e scattering time which is consistent with the 
spin Coulomb drag^£ observed in the Ref. @. Therefore, 
in the following discussion at finite T > 30 K, we can 
safely neglect the impurity scattering in our theory even 



though it can be trivially in corporate when extending 
the results to lower temperatures. 



B. Thermal average quantum kinetic equation 

In the equilibrium state, the Keldysh Green's function 
satisfies 2 ^ 

G*(E,k) = (G R - G A ) tanh(^-^). (19) 

When the quasiparticle approximation is valid, the 
Keldysh Green's function in the E — k space is still a 
peak even in the non-equilibrium state and has the form 

G K {E,k;T,R) = -27ri<5(E- e k )h, (20) 

where /i(k, R, T) is the distribution function, defined as 

dE 



h k {R,T) 



2ni 



Gu.e{^ t ) 



(21) 



In the linear response limit, the non-equilibrium distri- 
bution function takes the form 

h k (R,T) = -^-g(e,R,T), (22) 

where N is the density of state, 9 is the angle between 
k and the x axis, and 



g= / N de k h k (R,T) 



(23) 



is the thermal average distribution function. We also 
introduce the density operator 

d 2 k 

d0_ 
27 



MR, T ) 
J ^ f N h k (R,T)de k 
J ^g(9,R,T). (24) 



Multiplying by — 7p and integrating over efc and E on 
both sides of Eq. ©, the left hand side takes the form 



d T g + V fl ■ { - V, g} + i[b ■ a, g] + 



where 



(25) 



V = J -/'(e fc )y(k)de k = V(M), 
b = [ -f'(e k )b(k))de k = b(k, ¥, 9), (26) 
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and the right hand side takes the form 



»r . dE r 

N de k — L 



, dE dE 1 d 2 k' 
N de k — -^^—^ x 



'2n 2tt (2tt) 2 
5G K \D R - D A ){G R - G A ) 

fdE'd 2 k'^ K f dE R A 

x(^-^)(tanh(^) + eoth(^)) 
p(i?,T) x — . 



(27) 



Therefore, the kinetic equation in Eq. ((5J) is converted 
to the equation of the thermal average distribution 
g(8,R,T) and the density function p(R, T) which takes 
the form 



9t5 + Vr-^V.s} 



• i(5.*,$] + J- = ffi^2(28) 

7"ee ^~ee 

Eq. ((28]) gives the coupled kinetic equation of the ther- 
mal average distribution function g and the density ma- 
trix /5(R, T). However, the spin-charge dynamic equation 
normally is the equation of the density p. Therefore it 
is necessary to convert Eq. (|28[) to an equation only con- 
taining p. This is not easily done in the presence of SOCs. 
Mishchenko et o/.— treated the gradient term VG K , as a 
perturbation and developed a self consistent method to 
derive the diffusion equation to any order of V ' g in prin- 
ciple. However this method is restricted to the Rashba 
SOC and becomes harder when calculating the higher 
order gradient terms or considering the linear and cubic 
Dresselhaus SOC. 

Bernevig et aL— pointed out that in the ballistic 
regime ql > 1, where q is the spin-charge wave length 
and I is the mean free path, we have to consider the 
higher order gradient terms VG . This means that in 
the ballistic regime, the spin-charge dynamic equation is 
not dominated by the second order of the spacial differ- 
ential operators V 2 but one needs to consider the infinite 
summation over the gradient expansion. However, this 
can only be done for the very special case when a = /3i 
and /?3 = 0J£ For a generic SOC, it is a daunting task to 
evaluate the infinite summation. Therefore in this work, 
we abandon the idea of the gradient expansion of VG K 
and provide a different way to obtain the spin-charge dy- 
namic equation for the general SOCs. 



The third term on the left hand side of Eq. ([28]) has the 
form 

\bxa x + b y (jy,g} = 2i ((b x g y - b y g x )a z - b x g z a y + b y g z a x ) . 

(30) 

This couples the different spin components and generates 
spin precession. The second term on the left hand side 
of Eq. (|28|) contains the SOC velocity operators v so — 
dH so /dk which gives 

{\\<j y d x + \2<Txd y ,g} 
= (Xid x g y + \2d y g x )<r + \\d x g c cr y + X 2 d y g c a x .(31) 

This couples the charge and spin through the finite Vg 
which indicates non-uniform charge or spin distribution 
in real space. The other terms of Eq. (|3U|) do not cou- 
ple spin or charge components. If we multiply oi where 
% = 0, x, y, z and calculate the trace, using the fact that 
Tr((Ti(Tj)/2 = Sij, Eq. (|28]) can be rewritten as 

(&r + - + — ■ V)g c + \\d x g y + \2d y g x ) er = — ctq, 
t m J t 

(d T + - + — ■ V)g x + \ 2 d y g c + 2ib y g z J a x = —a x , 
t m J t 

1 fc - \ P 

ipT H 1 ■ V)g a + Xid x g c - 2ib x g z a y = —a y , 

t m J t 

Ik - - \ p z 

(d T H 1 • V)g z + 2i(b x g y - b y g x ) a z = —a z , 

t m J t 

(32) 

where k is the thermal average of the momentum. Above 
we have multiplied again by the corresponding matrix to 
remind ourselves of which component belongs to which. 
Hence, we can obtain the 4x4 kinetic equation of the 
coefficients g c i x>y>z ) and p c i x>y>z ) which takes the form 



K 



( 9c 

gx 

9v 
V 9z 




(33) 



where 



K = 



i\\q x T 




-i\2q y r i\\q x T 
Q 

o n 

2b y T —2b x r 



~2b y T 
2b x t 



(34) 



C. Spin-charge density dynamic equation valid in 
both weak and strong SOC regime 

The thermal average distribution function and density 
matrix can be generally written as 



9 z 
P : 



ScOo + g x <?x 

- pc&0 + Px&x 



g y o 

-PyC 



- PzOz- 



(29) 



Cl = 1 — iuT + q ■ vt, and 9 is angle between k and the 
x-axis. Here we have Fourier transformed 8t and d x r y -\ 
to — iuj and iq x ( y ) which are the frequency and wave- 
vector of the spin polarization wave in the experiments. 
To further obtain the spin dynamic equation, we simply 
multiply K^ 1 on both sides of Eq. (|2"g)l and integrate out 
the angle 8. The inverse of K is easily obtained for the 
generic SOCs because it is just the inverse of a 4 x 4 
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matrix. This provides a way to derive the spin dynamic 
equation in the presence of the general SOCs. 

Before discussing the spin dynamic equation, we would 
like to emphasize two advantages of our method. First, 
the inverse of the matrix K is equivalent to the infinite 
summation of the gradient expansion but much easier to 
be calculated exactly. In the Sec. IV. it is shown that 
in the case of a = /?i and /3 3 = 0, we obtain the same 
spin dynamics modes to those obtained by the infinite 
summation of the gradient expansion of VG K Second, 
we need not to guess the form of the non-equilibrium 
distribution function g before we obtain the spin dynam- 
ics equation. The non-equilibrium distribution function 
is often expanded on spherical harmonics (Eq. (7.79) in 
Refill) before we solve the kinetic equation. Normally in 
the system without SOCs, it is enough to stop the expan- 
sion on the first order of the spherical harmonics. How- 
ever, in the presence of the SOCs, especially consider- 
ing the cubic Dresselhaus SOC, the distribution function 
may contain the spherical harmonics up to third order. 
Therefore it is very difficult and complicated to guess the 
correct form of the non-equilibrium distribution function 
by expanding it to the right order in spherical harmonics. 
Our method does not have this difficulty and is approx- 
imate to considering all possible spherical harmonics in 
the non-equilibrium distribution function. 

To simplify our discussion of the spin dynamic equa- 
tion, we consider the spin wave vector only along the 
x-direction and take q y = 0. Because q <C fcp, the spin- 
charge coupling terms i\i(2)Qx(y) T are much smaller than 
the spin-spin coupling terms iA 1 ( 2 )fcr in Eq. and we 
neglect the spin-charge coupling and only focus on the 
spin-space of Eq. (l34l) . In this case, using the definition 
Eq. (|24|) the spin-charge dynamic equation of the density 
coefficient p c (x,y,z) can be obtained as 



(35) 





where D — f tt-K c 1 and K a 1 takes the form 



KZ 



I n 2 
\ 



Ablr 2 



Ab x b y r 2 
-2b v TVl 



4b x b y r 2 
Cl 2 + 4&Jr 2 
2b T rQ 



2b y rn 
-2b x rQ. 

n 2 



Ab 2 r 2 n 



(36) 



Eq. (|35[) describes the spin dynamics in the frequency- 
momentum space at finite temperature and any strength 
of the SOC. This is the key result of our theory. 



III. SPIN DYNAMICS IN THE WEAK SOC 
REGIME AT FINITE T 

In this section, we focus on the spin dynamics in the 
weak SOC regime which is found at temperatures above 
35 K in the experiments considered. 15 



A. Only DP spin relaxation mechanism 

In the regime where Ai(2)A;r <C 1, ql <C 1, defined as 
the weak SOC regime, the spin-charge dynamic equation 
can be written as 



1 



Px 




= 0,(37) 



where 



o 



^2 _ _ 

—i(a + 0)fl so q x -i(a-0)fl 
a 



i(a + 0)tt so q x 
0)1 
fit 





n 2 so (±-&0) i(a-0)fl so q v I , 



a 



= 



y/a 2 + Wi - ft) 2 + 0. 

01 - 03 

^/a 2 + (0 1 -0 3 ) 2 +0 : 



(38) 



Qx{y) = qx(y)VT, Q SO 



fla T and f2 s 



2^/ol 2 + (0i — 0s) 2 + 2 kp is the spin precession fre- 
quency due to the SOC. The first term on the left hand 
side of Eq. (f3"T| is the normal diffusion equation with- 
out SOCs. The second term corresponds to the average 
torque exerted by the SOC in the momentum space which 
gives the correction of the diffusion equation due to the 
generic SOCs. The detailed calculation of the elements 
of the matrix D is shown in App. [C] 

Our spin diffusion equation, Eq. (l37l). is equivalent to 
the spin diffusion equation in Ref. Il2| at zero temperature 
where k — kp. However only considering D-P mechanism 
is not enough to explain the temperature dependence of 
the enhanced-lifetime of the SHM^i 5 - which first increases 
with increasing temperature from 5 K to 75 K and then 
decreased with further increasing temperature. Because 
the density of states is a constant for the 2DEG, the 
chemical potential shift at finite temperature is of the or- 
der of (f^) 4 = 0.0039 when T = 100 K according to the 
Sommerfcld expansion. Our numerical calculation gives 
Ef = 398.15 K at T = 100 K which is consistent to the 
analytical prediction. Assuming that x(E) = —f'(E), at 
T = 100 K we have 



k= / dE- f'(E)k = 0.98k F , 



k 3 



dE - f(E)k 3 



l.lk 3 F , 



(39) 



which proves that the thermal average k and k 3 are not 
changed too much from K to 100 K. Therefore the in- 
crease of the thermal average cubic Dresselhaus SOC is 
not sufficiently strong to account for the non-monotonic 
T-dependence of the enhanced-lifetime of the SHM. An- 
other evidence for this statement is from Fig. 3 (c) in 
Ref. [H. The mobility is reduced to avoid the ballistic 
crossover and the spin life time is measured in five dif- 
ferent temperatures as a function of spin polarized wave 
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vector q. It is found that at q = 1.26 x 10 4 cm -1 , which 
is close to the SU(2) points, the spin life time is the max- 
imum at the lowest temperature 5 K and minimum at 
the highest temperature 250 K. When q is away from 
the SU(2) point such as q = or q = 2.5 x 10 4 cm -1 , the 
spin life time is the minimum at the lowest temperature 
5 K and maximum at the highest temperature 250 K. As 
a result, when q is far away from the SU(2) point, the 
enhanced-lifetime of spin helix mode matches the descrip- 
tion of D-P mechanism that the spin life time is inversely 
proportional to the momentum scattering time. When q 
is close to the SU(2) point, the turning point, at which 
the spin life time changes from increasing to decreasing 
with increasing temperature, is lower than 5 K. The 
turning point at such low temperature can not be the 
consequence of the thermal average of the cubic Dressel- 
haus SOC.-^ 



form 

V imp (r e ff) = V imp (r) - j(V r V imp (r) x k) ■ er, (41) 
and the Coulomb interaction now has the form 

V e - e (r le f f ,r 2 eff) = V e - e (r 1 ,r 2 ) ~ j{V ri V e - e x ki) ■ er 

- 7 (V r2 K-e x fe 2 ) ■ CT, (42) 

where the subscripts 1 and 2 of the space coordinates rep- 
resent the two interacting electrons respectively. There- 
fore the spin-orbit coupled scattering potential has the 
form 



V so (r) = -7(V P V(r) x k) ■ a, 



(43) 



where V(r) = (Vi mp + V e - e ) is the spin independent mo- 
mentum scattering potential. The spin relaxation rate 
due to the E-Y mechanism in the 3D bulk material has 
the form 31 



1 



B. EY mechanism in the 2DEG with the e-e 
interaction 

Because D-P mechanism is not sufficient to account 
for the T-dependence of the SHM life time, we have to 
consider the E-Y mechanism which is proportional to the 
momentum scattering time, opposite to the case of D-P 
mechanism. There are two processes involved: that of 
Elliott 28 and that of Yafet 29 . In the Elliott process, the 
scattering potential is spin independent. The spin flip is 
due to the SOC on the Bloch state, say the admixture 
of the Pauli up and down spins through the coupling be- 
tween the conduction and valence bands. In the Yafet 
process spin flips come directly from the scattering po- 
tential with the well known form 4m ^ c ^ (VV x P) ■ er so 
that the potential alone couples opposite spins. In the 
GaAs/AlGaAs 2DEG, as shown in App. El the Elliott 
mechanism is much larger than the Yafet mechanism. 
As a result, in the following discussion, we only consider 
the Elliott processes. The Elliott mechanism can be de- 
rived from the unitary transformation matrix based on 
the Lowdin partitioning^ The coordinate r after uni- 
tary transformation takes the form 3 ^ 



P 2 . 1 
r eff =r-—{ 



— -rr)k x er = r — 7fc x er, 

3 y E 2 (E a + A a ) 2 ' 1 



P 2 , 1 



3 y E 2 (E + A ) 2 



)k x er . 



(40) 



TEY (efc) 



Ar 



E a + A c 



Ea 



1 



(44) 



where the numerical factor A is of the order of 1 and 
dependent on the scattering mechanism™ such as e-e 
interaction 3 - 2 - and T(e^) is the momentum scattering time 
at energy e&. Because the non-equilibrium electrons dis- 
tribute around the Fermi surface within the energy range 
fcflT, the e-e scattering time is independent on the en- 
ergy efe according to Eq. (fT5)) and is labeled as r. In the 
2DEG, the momentum k z is quantized. In this case, al- 
though the average of (k z ) — 0, (k 2 ) w (^) 2 which gives 
the linear Dresselhaus SOC For the same reason, when 
considering E-Y mechanism, we also assume (k 2 ) ~ (^) 2 
where d is the width of the quantum well. Therefore, the 
E-Y mechanism in the 2DEG can be written as 



1 



1 



TEY,x T EY,y 
1 



A(l + 
4(fc 2 



TEY,-, 



A 



2(^ 2 ) , ( A 

An 



k\ \E + A a 



2 / > 2 , 

efc \ 1 

En 



kl VSo + A 



2 / x 2 -, 



E 



, (45) 



where TEY.x(y,z) are the spin relaxation time for the spin 
polarization along x, y, z direction respectively due to the 
E-Y mechanism. Appendix [A] gives the detailed deriva- 
tion of the E-Y mechanism in the systems considered 
here. 



Temperature dependence of the spin relaxation 
modes 



where P = ^- < S*|V|i? >, \S) is the s-wave like local 
orbital state, A — X,Y, Z are the p-wave like local orbital 
states, Eq is the band gap between T§ conduction band 
and r^" valence band and An is the SOC gap between 
the bands of and Tg i^l Any coordinates dependent 
potential will be modified by the second term in Eq. (|40| . 
For example, the impurity scattering potential has the 



In this sub-section, we show that the temperature de- 
pendence of the spin relaxation modes, especially the 
enhanced-lifetime of the SHM, is the consequence of the 
competition between the D-P and E-Y mechanisms. By 
adding the E-Y mechanism on the Eq. (|3"7| . the spin- 
charge dynamics containing both D-P and E-Y mecha- 
nism takes the form 



8 



-2 



f n x } 






*( 







nfod + fi /8) + «i i{& + p)h so q x 



—i(a + f3)Sl so q x -i(a - /3)fl so q v 









) : 


= 0. 


(46) 









where 



Kj 1 =A 



4(fc z 2 ; 

kp 



A 



So + Ao 

2 



So 







S + A 



So 



(47) 



Here, according to Eq. (|3T)1) . we take the thermal average 
of the momentum k = kp. To simplify our discussion, we 
assume that the spin polarization is uniform in y direc- 
tion and non- uniform in x-direction with finite q x . In this 
case the spin polarization along y-direction is not coupled 
to other two components of spin polarization. Based on 
the Eq. (|4"6")l . the eigen-mode of the spin polarization 
along y-direction, has the form 



(48) 
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FIG. 4. The reduced-lifetime and enhanced-lifetime of the 
SHMs. The x axis is the spin wave vector normalized by 
Q = 4m e //a and the y axis is the spin life time normalized 
by the spin precession period T s — Q so t = 0.1 is taken 

to satisfy the weak SOC condition and fia/fii = 0.16 from 
the experimental data^ The maximal enhanced-lifetime is 
still very close to Q although the cubic Dresselhaus SOC is 



where uj y is the normalized eigen-frequency of spin po- 
larization along y-direction. At the SU (2) symmetric 
pointjii where a — (3\ and (3 3 = 0, (1/2 — a/3) = and 
the effect of the D-P mechanism on the spin polarization 
along y-direction is zero. This is consistent to the fact 
that the effective magnetic field due to SOC in this case 
is always along the y-direction. 

Next we focus on the other two spin dynamic modes 
which have the form 



± 



\ 



n? (f -\) 



— and k\ = — 



(49) 



where Co± are the normalized eigen frequencies of the 
reduced- and enhanced-lifetime of the two SHMs^ re- 
spectively. The maximum enhanced-lifetime of the SHM 
happens when a = /3%, fa = and q = 4m e ffa = Q. 
In this case, we have a — & — n 1 - /™ 2 -1- # 2 — xl 
q = ArrieffavfT = \/2Q 



Y and 



,r = \/2Q so . However, in the 
material such as GaAs, the cubic Dresselhaus SOC is the 
consequence of bulk inversion asymmetry (BIA) and is in- 
evitable. When considering the cubic Dresselhaus SOC, 
the maximus spin life time is observed at /?i — /?3 = Q- 15 
which is consistent to the theoretical prediction at zero 
temperature*^ When 03 <C a, the maximum enhanced- 
lifetime of the SHM is still at Q = Am e ffOt shown in 
Fig. I As a result, the Rashba SOC strength can be 



detected by the relation a = Q /4m e f / which gives 



2 v /W(/3 1 -/3 3 ) 2 +/3 3 2 ^ 
2\f2ak F = 



Qk F 



V2r, 



(50) 



■eff 



Taking q x — 4m e ffa and substituting the above rela- 
tion to Eq. (l48l) . we obtain the reduced- and enhanced- 
lifetime of the SHMs as 



iuj± = — 

T 



ni 



K 2 T, 



(51) 



7 a 2 f ,a 2 1, «? 



8<5 2 



(52) 



where ± is corresponding to the reduced- and enhanced- 
lifetime of the SHMs respectively. 

We plot the enhanced-lifetime of the spin helix mode 
in Fig. [1] and compare our theoretical calculation with 
Fig. 3 (a) in the Ref. [TH The system we are consider- 
ing is the GaAs/AlGaAs quantum well in Ref. O Its 
width is llnm which gives hlh = {k 2 F )/^{k 2 z ) = 0.16^ 
Therefore, in this system, j3i — P3 — a 3> As and Q so is 
estimated to be 0.356 THz based on the Eq. §5U§. The 
band gap E = 1.519 eV, A = 0.341 eV. To fit the 
experimental data, Fig. [TJ the only fitting parameter we 
choose is A — 4.0 which is the order of 1 and consistent 
to the estimate in the Ref. [3l|. 

Based on the Eq. (f5"Tj) , both the reduced- and the 
enhanced-lifetime of the SHMs have a maximal spin life 
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time as a result of the competition between the D-P 
and E-Y mechanism. The critical momentum scatter- 
ing times rf- when the spin life time reach its maximum 
take the value 



-i- = 0.029ps, 

K+ 




-i = 0.51ps. 



(53) 



This is the result of the competition of E-Y and D-P 
mechanisms. It is noted that r + <C t_ which indicates 
the reduced-lifetime needs much higher temperature than 
the enhanced-lifetime to reach its maximum. Actually 
t+ = 0.029 ps is smaller than the minimum momentum 
scattering time 0.04 ps estimated from the spin life time 
of the enhanced SHM is 100 ps at T = 300 K based on the 
Eq. ([52"|). This is why the reduced-lifetime of the SHM 
increase monotonically with increasing T. 



where il so = 2\ikpT. The eigenmodes satisfy the equa- 
tion 



(1 - D u ) 2 + Dj 3 



(l-iw) 2 + (q x -n so ) 2 



-1 + y (1 - iCjf + (q x + n so j 2 

X , : = 0, 



(1 - iCo) 2 + [q x + fl 



which gives the spin relaxation eigenmodes as 



(55) 



<w = i±^i-(5 B ±n„) 2 , (56) 

and is consistent to the result in Ref. fl3l 

B. Only Rashba or linear Dresselhaus SOC 



IV. SPIN RELAXATION IN THE STRONG SOC 
REGIME 



In this section we show that our theory can reproduce 
the theoretical results of the spin relaxation in the case 
of a = Pi, = and in the case of only Rashba or 
linear Dresselhaus SOC in the strong SOC regime at zero 
temperature i 13 ' 14 



The persistent SHM occurs in the case of a = /3i, 
P3 = and q is along x-direction. In this case, the el- 
ements of the matrix D, corresponding to the reduced- 
and enhanced-lifetime of the SHMs, in the spin-charge 
dynamic Eq. (|35l) take the form 



1 



£>n = D 



33 



2^/(1 -^) 2 + (& + O so )2 
1 



2^/(1 - iu) 2 + {q x ~ n so ) 2 
i 



Dis = -Dsi = 



2^/(1 - lQ) 2 + (q x + n so ) 2 
i 



(54) 



2 J (1 - iu) 2 + {q x - fl so ) 2 



In the presence of Rashba or linear Dresselhaus SOC, 
we only discuss spin relaxation of the uniform spin po- 
larization because, as far as we are aware of, there is 
no analytical form of the spin relaxation eigenmodes for 
the non-uniform spin polarization in the presence of only 
Rashba or linear Dresselhaus SOC. For the uniform spin 
polarization, the off diagonal elements of the dynamic 
matrix D in Eq. (|35p are zero because the angle average 
of the effective magnetic field b is zero. Therefore, we 
only need to focus on the diagonal terms which have the 
form 



Du = D 2 



fr + n 2 so /2 



D^ = 



w + n 2 



(57) 



The spin relaxation mode of the spin polarization along 
the z-direction satisfies the equation 1 — D33 = and is 
solved to have the form 



1 



1 



1 - 4«L , 



(58) 



which is identical to our previous theoretical result^ in 
the zero temperature limit which has been shown to be 
consistent with the experimental observation. 17 



V. CONCLUSION 

We have developed a consistent microscopic approach 
to explore the spin dynamics in the presence of SOC of 
arbitrary strength, impurity scattering and e-e interac- 
tion at finite temperature. The e-e interactions are in- 
troduced by using the non-equilibrium Green's function 
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formalism. To create a full understanding of the spin- 
dynamics we incorporate both the DP and EY mecha- 
nism. We have shown that near the SU(2) symmetry 
point, because the D-P mechanism is suppressed for the 
SOC enhanced SHM, the E-Y mechanism dominates in 
the high temperature for this mode. By choosing the 
reasonable parameter of the E-Y mechanism, we show 
that our theory of enhanced-lifetime of the SHM matches 
the experimental data quantitatively. Our theory is also 



shown to be able to recover the several previous theoret- 
ical results, which are only valid in the non-interacting 
system and most focus on the zero temperature limit. 

We acknowledge Chia-Ren Hu, Artem Abanov, 
Alexander Finkel'stein and Valery Pokrovsky for very 
helpful discussion and support from DMR-1105512, 
lonr-n000141TH 0780. NHARP, and NSF-MRSEC DMR- 
0820414. 



Appendix A: Elliott- Yafet mechanism 

In this appendix, we discuss more extensively the spin relaxation due to the Elliott- Yafe t 28 > 29 mechanism in the III- 
V semiconductor quantum well. There are two processes involved: the Elliott process a.ndM and the Yafet process.—. 
In the Yafet process, the spin flip is due to the intrinsic spin-orbit coupling of the scattering potential which has the 
well know form 



HI 



Arrive 2 



a ■ {k x W), 



(Al) 



where mo is the electron mass in the free space and V is the scattering potential which can be impurity, e-ph and 
e-e scattering. Here to estimate the strength of the Yafet process, we assume V(r) is a spherical potential which is 
independent on the direction. In this case, the SOC Hamiltonian is Eq. (|A1|) is simplified to be 



HI 



Am 2 c 2 



cr{k x e r d r V(r) = — 



d r V{r) 



AttIqC 2 



(r x k) ■ cr = — 



d r V(r) 



AninC 2 



h 2 d r V(r)j(j + l)-l(l + l)- S (s + l) 
Arrive 2 r 2 



(A2) 



where L is the orbital angular momentum, j(l,s) is the quantum number of the total angular momentum (orbital 
angular momentum, spin). In the of the III-V semiconductors, j = 3/2, I = 1 and s = 1/2. Therefore the 

.d r V(r).l h 



expectation value of H™ trt in T 7 can be written as 



h v 1 = (j = 



,l = l\H so \j = -,l = l) = - 



AiuqC 



2„2 



Arrive 



2„2 



V 



Am 2 c 2 (a /2) 2 2 

T7 r ? 



(A3) 



tV and 



where do is the lattice constant. Here we assume that (1/r 2 ) = 4/a 2 ,. Similar, in the TU band, H S 7 Q = 

m u- "U 

in the Tg, Tf s l = 0. 

In the Elliott process, the spin flip is due to the spin-orbit interaction (SOC) on the Bloch state, say the eigenstate 
of the electron in the conduction band is not the spin eigenstate. However the momentum scattering potential is V 
which is spin independent. From this viewpoint, the Elliott process has the same origin as the Rashba SOC in the 
III-V semiconductor quantum well. 

The difference lies the fact that in the Rashba SOC, the electric field along z-direction break the inversion symmetry 
and the first non-zero term of Rashb SOC is linear dependent on the wave vector k of the conduction electrons which 
make the spin-up and spin-down states with the same momentum k non-degenerate. In the Elliott processes, the 
scattering potential does not break the inversion symmetry, its nonzero terms are proportional to k 2 and spin-up and 
spin-down electrons with the same momentum are still degenerate. 

Based on the 8x8 Kane model, the effective Hamiltonian of the electrons in the conduction band of III-V semi- 
conductors such as GaAs has the form 



H n = — 



P 2 
3 

El 

3 



2 
1 



1 



£o + A 
1 



p2 

k 2 - — 

3 

h 2 



l 

^0 



1 



(E + A ) 



(Eo + Ao) 2 
er • W x k, 



cr • W X k 



(A4) 



'o l 



'0 



where V is the spin independent scattering potential, P = ^ < S\W\R >, \S) is the s-wave like local orbital state, 
R = X,Y, Z are the p-wave like local orbital states, Eq is the energy gap between Tg" and T$ bands and Aq is the 
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energy gap between and bands^I. The second term in Eq. (IA4[) is corresponding to the Elliott spin relaxation 
mechanism and vanish when the SOC gap A = 0. The third term in Eq. (|A4[) is corresponding to the Yafet spin 
relaxation mechanism which will not vanish when the SOC gap Ao = 0. These are consistent to the characters of 
Elliott and Yafet mechanisms. For the Elliott spin relaxation mechanism, the scattering potential is spin independent, 
therefore only the SOC gap Ao can provide the necessary SOC to relax spin. In the Yafet mechanism, the scattering 
potential itself contains SOC which can relax the spin. Therefore the SOC gap Ao is not necessary in this case. 

Now, let us compare the strength of the Elliott and Yafet mechanism in the GaAs 2DQW, where P — 10.493 eVA, 
and the band gaps Eq = 1.519 eV, Ao = 0.341 eV^£ Therefore, in the III-V semiconductors, the ration between the 
Elliot and Yafet mechanism is of the order of 





1 1 








1 




K 2 


(E +A ) 2 _ 


2m' 2 j c 2 a' 2 j 



= 4.84 x 10 5 , 



(A5) 



which indicates that the Elliot spin relaxation is much larger than the Yafet spin relaxation in a III-V semiconductor. 
Therefore, in the following derivation, we only focus on the Elliot mechanism. The first term in the Eq. (IA4|) is 
the normal kinetic energy and the second term is the SOC terms In 2D quantum well, by applying an electric 
potential along z direction V z — —eEz, the inversion symmetry is broken and (V14) is nonzero which gives the 
Rashba SOC. However, in the Elliott process, the potential V is from the impurity scattering, clectron-phonon 28 and 
electron-electron interaction 3 - 3 ., the average, (VT^), is zero. Therefore the first nonzero term of the Elliott process is 
proportional to ((W)) 2 . The SOC potential can be written in the momentum space as 34 



V so (k,k') 



ijV(k x k') 

hv{{k x k' y 



k y k' x )a z + (k y k' z - k z k' y )a x + (k z k' x - k x k' z )a y } = i{V z a z + V x a x + V y a y ), (A6) 

) 



where k' and k are the momentums of the electrons before and after scattering respectively, 7 = -^-(-gs- — (E +A y A 
and V is the momentum scattering potential. 

Before we calculate the Elliott mechanism, let us connect the quantum kinetic equation to the continuous equation 
by integrating the momentum k. To simplify our argument, we consider the case where the SOC due to the inversion 
asymmetry is zero. Eq. © after the integral has the form 

d 2 k 



d T p{E, R,T) + V- J(E, R,T)= [ \{^ R G K + ^ K G A ) - (G R T, K + G K ^ A )} 

J {2tt) z 



(A7) 



Eq. (|A7I) is the continuous equation in the spin-i basis. Multiplying by a m /2 on both sides of Eq. (|A7|) and taking 
the trace over the spin space we have the traditional continuous equation 

d TPm + V • J m = j ^Tr^p [(£ R G K + S K G A ) - (G R £ K + G K E A )] , (A8) 

where m — 0, x, y, z are corresponding to charge, spin-cc, sp'm-y and spin- z respectively. If the charge or spin is 
conserved, the scattering term on the right hand side of Eq. (|A8|) is zero. Now let us substitute the SOC potential 
into the scattering term. There are four terms in the collision integral and we first focus on the two terms containing 
the retarded Green's function or self energy which has the form 

d 2 k' a n 



*k_ 1 

(2tt) 2 2 



S R G K - G K S 



RvK 



>} 



Tr 



= Tr 



Tr 



(27T) 2 
d 2 k 



2 

d 2 k' 



(2tt) 2 (2tt) 2 
d 2 k d 2 k' 



V kk ,G R (k')V k , k G K (k) - G R (k)V kk ,G K (k')V k <k 
(^V kk , - Vkk ,^)G R (k')(V k , k )G K (k) 

(G R (k')V k ' k G K (k)Vkk< - G R (k)V kk ,G K (k')V k , k )(^9) 



(2tt) 2 (2tt) 2 2 

The last term on the right hand side of Eq. (IA9|) is zero because the symmetry of k and k' in this term. The first 
term on the right hand side of Eq 

d 2 k' - 



9l can be written further as 



Tr 



Tr 



Tr 
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d 2 k ; 

d 2 k' 



(^V kk , - Vkk ^)G R (k')(V k , k )G K (k) 



[^,VV]G R (k')(V k , k )G K (k)^ =Tr 
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-G R (k')G K (k) 



Tr • 



d 2 k' 



4(> ;vL,vt 



kk' v k'k^ 



-G R (k')G K (k) KA10) 
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For the charge dynamics, we need to take m = 0. Because cto commutes with any scattering potential operator, 
Eq. flMOj) is always zero which indicates that the charge is always conserved. For the spin dynamics, m — x,y,z. 
Generally speaking, c x (y >z ) does not commutes with a spin dependent scattering potential and the Eq. (|A10|) is nonzero 
which indicates that the spin will decay. Note that in the 2DEG, k z is quantized. Therefore, taking m — x,y,z we 
have 



Tr / 4^4(5: vi k ,Vi, k )fG«(kOG K (k) = -^(1 + ^)Gf = + ) 

j#x J P F F 



Tr 



Tr 



(2 , )2 4(EvLvk)^ ( ko G K ( k )| =-^ ( i + f) G f = -^(i, fc| , Uo + Ac 

d2k/ lf y v j V j )^G R (k')G K fk)l- 4 l!M:lM G ^- 4fc,V A \ 2 / £fc \ 2 Gf 





\ 2 Gf 












I 2t/ 




(All) 



where fc z — n/d and ci is the width of the quantum well and A is the order of unity. The other two terms, containing 
the advance Green's function in Eq. (|A8[> . gives the same result to the Eq. (|A11[) . The average, (k z ), is nonzero 
which is also the reason why there is a linear Dresselhaus term in the III-V semiconductor quantum well. Therefore, 
although average (k z ) is zero, the components of Elliott mechanism which are proportional to (k z ) are still finite. As 
a result, the Elliott spin relaxation rates have the form 

i _ i =,4(i+ W A ° w*Vi 



tey.x t ey . v k% \E + A J \E 

J-=A^(*Y(«lY±. = !±. (A12) 

TEY,z kp \Eq+A J \E J Tp Tp 

The same result can be obtained from the traditional definition of the spin decay rate due to the admixture of the 
Pauli spin-up and spin-down in the eigenstates of the conduction electron. From the 8x8 Kane mode we can obtain 
these eigenstates as 

-IP, ,3 3, l~2 P .3 1. 1 P , .3 1. 

^ - + T2E~o k+l r 2 } + V alb** 1 * 2> + TeE-o^r-J 
-i p , .i i. -i p , ,i i. 



V3^o + A Z| 2'2 / V3^o + A '2' 2 7 

-IP. ,3 1, [2 P , .3 1. 1 P , .3 
^ = |5 ' ;> + 7!^ fc+l 2' 2> + V 3^2' ~2> + Tl^- 1 ^ 
_3 -1 P 1 1 1 P 1 1 

2 ) + V3i?o + A fc+l 2'2 ) + V3^o + Ao^'2' 2 h 



The transition amplitude between these two states with different wave vector k are proportional to 
(*cfc'jnr)|¥<*, t ) - (-^( fc - fc - - fc - fc i) + 5 {Eo + A ) 2 (fc ^ L " ^ 

/ 1 p2 p2 

= V Q [--(k z k'_-k_k' z )(- 



(A13) 



3 V "*'"" "—'"Eg (£o + A )2' 
= -V 7(fc.#- - fc_*£). (A14) 

Therefore the spin decay rate is proportional to 

— oc V^ 2 {k z k'_ - k-k' z )(k z k' + - fc+A^) cx 7 2 fc/ — , (A15) 

where l/r p is proportional to the Vq. Eq. (IA15[) is equivalent to the Elliott spin relaxation time we obtain from the 
projected Hamiltonian into the conduction band. 
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Appendix B: How the e-e interaction vanishes in calculating the charge conductivity 

When deriving the quantum kinetic equation, the collision integral contains many kinds of momentum scattering 
channels; one dominant at finite temperature is the e-e interaction. However as we know that when calculating the 
conductivity, the e-e interaction should vanish in this case. In this appendix, we show that how the e-e interaction 
naturally vanish in calculating the conductivity and, therefore, they only affect the spin relaxation channles. To 
simplify our discussion, we start from the classical Boltzmann equation without SOC 

drfk, + ^{V fe , -VrAJ - eE • V fe / fel = J U kuK {f kl - f k{ ) 
+ ( (Tfciki, tf^Wto, fcj, fcs, h) [/(fc 1 )/(fc 2 )(l - f(k[))(l - f(k' 2 )) - f(k[)f{k' 2 ){l - f(h))(l - f(k 2 ))} , 



(Bl) 

where n is the dimension of the system, /(fc) is the Fermi distribution function, U kl _ k i i is the impurity scattering rate 
and W(ki, k[; k 2 , k' 2 ) is the e-e scattering rate. If we exchange /i, f 2 and /{, / 2 , the right hand side of Eq. (|B1I) is 
unchanged. This means the electron 2 satisfies the same Boltzmann equation as electron 1: 

d T fk 2 + i{V fc , -VR/feJ - eE • V k f k2 = J V(f k2 - f K ) 

d n k[d n k 2 d n k' 2 W(k u k 2 , k 3 , k 4 ) [/(fc 1 )/(fc 2 )(l - f(k[))(l - f(k' 2 )) - f(k[)f(k 2 )(l - /(fcx))(l - f(k 2 ))} , 

(B2) 

The charge current operator is defined as 

J = J d n k, (-e^r) f(ki) = J d n k 2 (-e^r) /(fe), (B3) 

where m* is the effective electron mass. Since we are interested in the dc-conductivity, the system is uniform and 
independent on time. Therefore the first two terms on the left side of Boltzmann equation are zero. To get the charge 
current equation from the Boltzmann equation, we multiply the charge current operator on the both sides of Eq. (|B1[) 
and Eq. (|B2|) which gives 

f dkl ^l E -d k f 1 = [ dk!dhi^U(kuK)(fi- fi) 

J m J m 

+ J dkxd^dk^dk'^-WikxMXX) l/i/2(l - f[)(l - f 2 ) - f[f 2 (\ - /0(1 - / 2 )) . (B4) 



f dk 2 ^E ■ d k j 2 - / dk 2 dk' 2 -^U(k 2 ,k' 2 )(f 2 - f 2 ) 
J m J ui 



dhdk.dk'.dk'^Wikuk,, k[, k' 2 ) (A/ 2 (l - /0(1 - f 2 ) - f[f 2 (\ - h)(l - f 3 )) . (B5) 
m 

The left hand side of Eq. jB4} and Eq. (|B5|) have the form 

f e 2 ki( 2 \ f ,e 2 ki( 2 ) s f , e 2 ki( 2 \ x ne 2 

J dk 12 ^^E ■ a k f m = J dk m d k ■ (—m Efm) _ J dfcl(2)( a fe _M )jB/l(2) = _ 

(B6) 

where n is the density of the electron and j = J dk^f(k) is the electric current. The first term on the right hand 
side of Eq. (|B4I) and Eq. (IB5|) correspond to the impurity scattering and equal to — where Ti mp is the momentum 
scattering time due to the impurity. The second term on the right hand side is not easy to calculate because the 
momentums ki, k 2 , k[ and k 2 are not independent but correlated by the fact that the e-e interaction conserve the 
net momentum, i.e. ki + k 2 = ki + k 2 . However, if we calculate \{Eq. \BA\i + Eq <\B5\i ), the e-e scattering can be 
written as 

dfc 1 dfc 2 dfcid^ e(fcl + fca) W(fe 1 ,fe a , h)', k' 2 ) (A/ 2 (l - f[)(l - f 2 ) - fim - A)(l - f 2 )) = 
Zm 

dkMKw^k,-, fc 1; fc 2 )( e(fcl 2 + fc2) /i/ 2 (i - f[)(i - f 2 ) - ^y^/^a - A)(i - h)) 

0. (B7) 
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Therefore the charge current equation, \{Eq. (ji?4[) + Eq. (|i?5p ). has the form 

ne 2 j 

(B8) 

where e-e scattering term exactly disappear in this current equation as long as the e-e interaction conserves the net 
momentum which is not sensitive to the form of the e-e scattering rate W{k\, k'^, ki, k' 2 ). However, when deriving the 
density matrix equation, we do not multiply by the current operator and therefore e-e scattering can not be canceled 
as we did in the current equation. This is reasonable because the density does not care about the net momentum 
which is the key to get the current equation. Therefore although the original Boltzmann equation of the distribution 
function contains the e-e interaction, the electron-electron scattering will disappear when we calculate the charge 
conductivity. Hence the e-e interaction will not affect the conductivity but can contribute the spin dynamics. 



Appendix C: The derivation of the matrix element of the spin dynamic equation 

The denominator of Eq. (l3"rJl) can be expand to 

1 



Cl 2 + 4k 2 T 2 cos 2 9X\ + Ak 2 r 2 sin 2 9X\ 
« 1 + 2iQ - 2i(q x cos 9 + q y sin 6) - 3ql cos 2 9 - 3g 2 cos 2 9 - Ak 2 r 2 cos 2 9X\ - 4fc 2 r 2 sin 2 9\\ (CI) 

where oj = ujt and q x {y) = Qx(y)VT. Substituting Eq. (|Cip to Eq. (|28|) we have 

d9 2X 1 krcos9 



D 13 = -D 31 



2ir n 2 + Ak 2 r 2 cos 2 9X 2 + Ak 2 r 2 sin 2 9\ 2 



/dO 
—Aikrq x {a + ft - 2ft cos26») cos 2 9 = -2i(a + ft - (3 3 )k T q x (C2) 
2tt 



D23 — —D32 — 



-2A?fcr sin0 



2?T f} 2 + 4fc 2 r 2 cos 2 dX j + 4fc 2 T 2 sin 2 Q} 2 



f rlf) 

/ — AikrqJPx - a x + ft cos 29) sin 2 9 = 2i(ft - ft - a)krq y (C3) 
J 2tt 



D 33 = 



d,9 



/ 



2tt n 2 + 4k 2 r 2 cos 2 9X\ + Ak 2 r 2 sin 2 9X 2 
d9 



— (1+iQ-qt cos 2 9 - qi sin 2 9 - Ak z r 2 cos 2 9Xf - 4k 2 r 2 sin 2 9X 2 2 ) 
2tt 



1 + iu,- \{ql + q 2 ) - 4(a 2 + (ft - ft) 2 + ft 2 )fcV = 1 + iCj - \{q 2 x + f y ) - n 2 so r 2 



(C4) 



f d9 n 2 + AX 2 k 2 r 2 sin 2 9 

~ J 2^ Cl{Ct 2 + 4k 2 r 2 cos 2 9X 2 + Ak 2 r 2 sin 2 9X 2 ) 

~D 33 + [ ^4A 2 fc 2 r 2 sin 2 9 = D 33 + (i (a - ft) 2 + (a ~ ft) ft + ft 2 )4fc^ 2 
J 2tt 2 

= l+id- \(ql + q 2 y ) - 2(a 2 + (ft - ft) 2 + ft 2 )fc 2 r 2 - a(ft - ft) 

= 1 + iu, - \{q 2 x + q 2 y ) - ir! 2 r 2 4fc 2 ? r 2 a(ft - ft) (C5) 
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£>22 = 



de n 2 + 4\?,k 2 T 2 sin 2 e 

2^n 2 + 4k 2 r 2 cos 2 9X 2 + Ak 2 T 2 sin 2 OX 2 . 

I ^(i +iCj-q 2 x cos 2 9 - q 2 sin 2 6 - Ak 2 r 2 cos 2 QX\) 
J 2n 

1 + i£j - \{q 2 x + <f y ) - 2(a 2 + (ft - ft) 2 + Pl)k 2 r 2 + a(ft - ft) 

1 + iw - i(g 2 + g 2 ) - ^ 2 D r 2 + 4fc 2 r 2 a(ft - ft) (C6) 



where tt so = 2V« 2 + (ft - ft) 2 +/3 2 fc. 
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